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Abstract. We show that the global fluctuations of spectra of GOE and GUE matri- 
ces and their principal submatrices executing Dyson's Brownian motion are Gaussian 
in the limit of large matrix dimensions. For nested submatrices one obtains a limiting 
three-dimensional generalized Gaussian process; its restrictions to two-dimensional 
sections that are monotone in matrix sizes and time moments coincide with the 
two-dimensional Gaussian Free Field with zero boundary conditions. The proof is 
by moment convergence, and it extends to more general Wigner matrices and their 
stochastic evolution. 



Introduction. The fact that the global spectral fluctuations of a GOE or a GUE 
random matrix evolving under Dyson's Brownian Motion, are asymptotically Gauss- 
ian is well-known, see §4.3.3 in [AGZ] and references therein, and also [S] for a gen- 
eral f3 analog. On the other hand, it was shown in [B] that the global fluctuations 
of spectra of various principal submatrices of a single GOE or GUE matrix are also 
Gaussian. The goal of this note is to put these two statements together. 

We prove the asymptotic Gaussian behavior for submatrices of a class of stochas- 
tically evolving Wigner random matrices that includes Dyson's Brownian Motion 
for GOE and GUE. The proof is by the method of moments, and the argument is 
slightly more general than the one presented in [AGZ] for a single Wigner matrix. 

We also compute the resulting covariance kernel explicitly. In the case of nesting 
submatrices, it represents a three-dimensional generalized Gaussian process, where 
one dimension comes from the position of the spectral variable, the second dimen- 
sion reflects the size of the submatrix, and the third dimension is the time variable. 
When restricted to the two-dimensional sections that are monotone in matrix size 
and time variables, it reproduces the two-dimensional Gaussian Free Field (GFF) 
with zero boundary conditions. 

In the case of GUE, the appearance of GFF on monotone sections could have 
been predicted from the determinantal sctructure of the correlation functions [FF] , 
[ANvMl], and from the analysis of [BF] that showed how such a structure leads 
to GFF covariances in the global asymptotic regime. However, the complete three- 
dimensional covariance structure seems to be unaccessible via that approach for 
example because the spectra of the full set of submatrices evolve in a non-Markovian 
way [ANvM2]. 
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Wigner matrices. Let {Zij(t)}jyi>i^i^K and {l^i(t)}i>i.teR be two families of 
independent identically distributed real- valued stochastic (not necessarily Markov) 
processes with zero mean such that for any fc > 1 

max(E|Zl2(^)|^E|yl(^)|'=) < oo. 

Set c{s,t) = ^KYi{s)Yi{t) and assume that 

c(s,t)>0, c{t,t) = l, 
EZi2(s)^i2W = c{s,t), V.Zl^{s)Zl^{t) = 2c{s,tf + 1. 



Note that by Cauchy's inequality c{s,t) < ^Jc{s,s)c{t,t) = 1. We say that a 
function c(s, t) is admissible if it arises in this way. 

One possibility for the above relations to be satisfied is to take all {2~iYi{t)} 
and {Zij{t)} to be independent standard Ornstein-Uhlenbeck processes on M; then 

c(s, t) — exp(— |s — We will refer to this possibility as to Gaussian specialization. 
Define a (real symmetric) time- dependent Wigner matrix X (t) = [X{i,j \ i)]ij>i 

by 

Ztjit), i<j, 



X{i,j\t)=Xij,i\t) 



Yi, i = j. 



An Hermitian variation of the same definiton is as follows: Let {Z^j} jy^yi now 
be complex-valued (i.i.d. mean zero) stochastic processes with the same uniform 
bound on all moments. Denote d{s, t) = EFi {s)Yi {t) and assume that 

d(s,t)>0, d{t,t) = l, 
EZi2(s)Zi2(i) EE 0, V.Z^2{s)Z^) = d{s,t), E|Zi2(s)nZi2(t)|2 = d{s,tf + L 

We will also say that a function d{s,t) is admissible if it arises in this way. 

There is also a Gaussian specialization that corresponds to {Yiit)} and {2^'RZij{t)}, 
{2^'^Zij{t)} being independent standard Ornstein-Uhlenbeck processes on M; in 
that case d{s,t) = exp(— |s — t\). 

Define an Hermitian time-dependent Wigner matrix X{t) = [X{i,j \ t)]i,j>i by 

Ztj{t), i<j, 
Yi, i=j- 



Xii,j\t) = Xij,z\t) 



Under the Gaussian specializations, the matrix stochastic processes defined above 

are called Dyson's Brownian motions. Traditionally one distinguishes the two cases 
by a parameter /3 that takes value 1 in the the real symmetric case and value 2 
in the Hermitian case. The random matrices arising at a single time moment are 
said to belong to the Gaussian Orthogonal Ensemble (GOE) in the P = 1 case, and 
Gaussian Unitary Ensemble (GUE) in the (3 = 2 ease. 

The height function. For any finite set S C {1,2,...} denote hy Xb the x 
submatrix of a matrix X formed by the intersections of the rows and columns of X 
marked by elements of B. 
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The height function H associated to a time-dependent Wigner matrix X is a 
random integer- valued function on M x ]R>i x M defined by 

H{x,y,t) = ^ ^ {the number of eigenvalues of X|i^2,...,[y]}(0 that are > x}. 

More generally, let A = {a,i}„>i be an arbitrary sequence of pairwise distinct 
natural numbers. Then we define the height function Ha via 



HA{x,y) 



{the number of eigenvalues of -'l^{ai,...,or that are > x}. 



The first definition corresponds to A = N. 

The convenience of the constant prcfactor /97r/2 will be evident shortly. 

A three-dimensional Gaussian field. Let c{s,t) be an admissible function as 
defined above. Set H = e C | 9^ > 0} and introduce a function 

C : (H X R) X (H X R) -)> K U {+00} 



via 



C{z,s]w,i) = 



c(s, t) min(|zp, — zw 



c(s, t) m\n{\z\^, — zw 



= < 



/ 

1 


In 


c(s, t)z 


- w 


27r 


c{s,t)z 


- w 


< 

1 


In 


c{s, t)w 


— z 


27r 

V 


c{s,t)w 


— 1 



\z\ < \w\ 



\z\ > w . 



It is easy to see that for any {s,t) with c{s,t) < 1, C( • , s; • , t) is a continuous 
function on H x H. Note also that if c(s, t) = 1 then 



C{z,s;w,t) = - ^In 



z — w 



z — w 



is the Green function for the Laplace operator on EI with Dirichlct boundary con- 
ditions. Viewed as a function in {z,w), it represents the covariance for the two- 
dimensional Gaussian Free Field on El with zero boundary conditions. 

Proposition 1. For any admissible function c{s,t) as above, there exists a gener- 
alized Gaussian process on H x K with the covaraince kernel C(z,s;w,t) as above. 
More exactly, for any finite family of test functions fm{z) S Co(H x R) the covari- 
ance matrix 



cov{fk,fi)= / / fk{z,s)fi{w,t)C{z,s;w,t)dzdzdsdwdwdt, k,l = l,...,M, 
JmJm 

is positive-definite. 

Denote the resulting generalized Gaussian process by Gc{s,t)- 
A proof of Propositon 1 will be given later. 
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Complex structure. Let A be a sequence of pairwise distinct integers. The height 
function Ha{x, y, t) (or H{x,y, t) = H^{x, y,t)) is naturally defined on MxR>i xR. 
Having the large parameter L, we would like to scale (x, y) ^ {L~2x, L~^y), which 
lands us in R X ]R>o x R. 

Wigner's semicircle law implies that for any t G R, with L » 1, x ~ is , y ~ L, 
after rescaling with overwhelming probability the eigenvalues (or, equivalently, the 
places of growth of the height function in x-direction) are concentrated in the 
domain 

{(x,!/) gRxR>o I -2^<x<2^]. 
Let us identify the interior of this domain with H via the map 



^■{x,y)^^ + i^y- (I) • 

Its inverse has the form 

n-\z) = {x{zU{z)) = {21Si{z)M''). 

Note that this map sends the boundary of the domain to the real line. 

Thanks to Vl we can now speak of the height function Ha as being defined on 
H x R; we will use the notation 

H^{z;t)=HA{L^x{z),Ly{z),t), zeU. 
Note that we have incorporated rescaling in this definition. 

Main result. Let X be a (real symmetric or Hermitian) time-dependent Wigner 
matrix. We argue that the centralized random height function 

H^\z;t) -EH^{z;t), zGU, tGR, 

viewed as distribution, converges as L ^ oo to the generalized Gaussian process 
ee(«,t) with c{s,t) = ^EYi{s)Yi{t). 

One needs to verify the convergence on a suitable set of test functions. The exact 
statement that we prove is the following. 

Theorem 2. Pick r e R, y > 0, and k G Z>o. Define a moment of the random 
height function by 

/ + 00 
x''{H{Lix,Ly,T) -EH{Lix,Ly,T))dx. 
-oo 

Then as L ^ oo, these moments converge, in the sense of finite dimensional dis- 
tributions, to the moments of Gc{s,t) defined as 

Mr,y,k= [ {x{z))''gc(s,t){z;T)^^^^^^dz. 



4 



Monotone sections as two-dimensional Gaussian Free Fields. Consider a 
time-dependent Wigner matrix and assume that the function c(s, = f Eil {s)Yi (t) 
is continuous and that it has the following monotonicity property: For any s e M, 
c(s, t) is strictly increasing in t G (—00, s] and it is strictly decreasing in f G [.s, +00). 
In other words, as time distance between matrices grows, the correlation decays. 
Further, assume that c(s, t) ^0 for any s,t gR. 

Let c6 : M — > ]R>o and -0 : R — > R be a continuous nonincreasing and a continuous 
nondecreasing functions, and assume that for at least one of these functions the 
monotonicity is strict. 

Our goal is to consider the joint fluctuations of spectra of matrices 

(1) ^{l,2,...,[L0(t)]}(VXi)), t€R, 

where L ^ 1 is a large parameter. By Wigner 's semicircle law, the spectrum of 
such a matrix scaled by is concentrated on [— 2-\/0(t), 2^/(j)(t)]. 

The two extreme cases arc 4>{t) = const (the size of the matrices is fixed and 
the time is moving) and ip(t) = const (the time moment is fixed and the size of the 
matrices is changing). 

Let us choose a reference time moment to gR and introduce a map 



S : {(a;, t) e M X M I -2^/^ <x< 2^/^} 



as 



E{x,t) 




t < to. 



[ c(V(0,V'(io)) 

The continuity and monotonicity assumptions on c, 0, and tjj are needed for S to be 
a bijection. Hence, its inverse is correctly defined, denote it as = {x{Q, t{()). 

We can now view the height function H for matrices (1) as a function on H via 

H^io = H{Li ■ x{o, L ■ mo), mo)) ■ 

Our main result implies that the centralized height function 

viewed as a distribution, converges as L — >■ cxd to the Gaussian Free Field on H in 
the sense of Theorem 2. 

Moments as traces. Let us rescale the variable x = L~iu in the definition of 
Mr,y,k and then integrate by parts. Since the derivative of the height function 

H{u, [Ly],t) in u is 

—H{u, [Ly],t) = -^P^Y.^{u- A,), 

S— 1 

where {^s}i<s<[Ly] are the eigenvalues of X^i we obtain 

\s—l s—1 



k + 1 

We can now reformulate the statement of Theorem 2 as follows. 
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Theorem 2'. Let X(t) be a time- dependent (real- symmetric or Hermitian) Wigner 
matrix with c{s, t) = ^EYi {s)Yi (t) . Let /ci , . . . , fc^ > 1 6e integers and yi,.. .,ym € 
M>o, ti, . . . , G M. Then the m-dimensional random vector 

(Tr(4r,...,[^,^j}(tp)) - ETr(xJr,...,[^,^jj(tp))^ 

converges (in distribution and with all m.om.ents) to the zero mean m-dimensional 
Gaussian random vector (^p)^i with covariance 



2kpkq 

^?PS9 — n 
pTT 



3;2>0 ?!w>0 



X — In 
2tt 



{x{z))''^-\x{w))''^-^ 

c{tp, tq) min (?yp, yq) - zw 



c{tp,tq)mm{yp,yq) - zw 



dx{z) dx{w) 
dz dw 



dzdw. 



More general submatrices. In the spirit of [B], we will actually prove a more 

general claim that involves arbitrary sequences of symmetric submatrices of the 
Wigner matrix that are sufficiently well-behaved. The exact statement is as follows. 

Theorem 2". LetX{t) be a time- dependent (real- symmetric or Hermitian) Wigner 
matrix with c{s,t) = ^EYi{s)Yi(t). Let fci, . . . , fcr„ > 1 6e integers, ti,.. .,tm G M, 
and let Bi, . . . , be subsets off^ dependent on the large parameter L so that there 
exist limits 



lim 



\B„ 



>0, 



L—t-oo Li 

Then the m-dimensional random vector 



lim 

-L— >oo 



\Bpf^Bq\ 



p,q = l,...,m. 



(2) 



L-^ Tr(X^aip))-ETr(X^ntp)) 



converges (in distribution and with all moments) to the zero mean m-dimensional 
Gaussian random variable (^p)^i with the covariance 



(3) IEe,^. = ^ 



{x{z))''^-\x{w))'' 



9!z>0 9tu>0 



X — In 
27r 



zw 



C{tp, tq)bpq - ZW 



dx{z) dx(w) 
dz dw 



dzdw. 



Theorem 2" can also be viewed as the moment convergence of the centralized 
height function Ha{x, y, t) to a limiting generalized Gaussian process but we do not 
give further details here. The static variant of this convergence is discussed in [B]. 
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Proof of Theorem 2". The argument closely follows that given in Section 2.1.7 
of [AGZ] in the case of one set Bj = B, and the proof of Theorem 2' in [B] in 
the static case. One proves the convergence of moments, which is sufficient to also 
claim the convergence in distribution for Gaussian limits. 

Any joint moment of the coordinates of (2) is written as a finite combination 
of contributions corresponding to suitably defined graphs that are in their turn 
associated to words. This reduction is explained in in Section 2.1.7 of [AGZ]. The 
key fact in the real-symmetric case is that averages of products of powers of matrix 
elements that involve at least one matrix element with exponent 1 vanish. The 
time-dependent analog of this fact is that averages of products of powers of matrix 
elements taken at different time moments that involve one matrix element with 
exponent 1 at only one time moment vanish. This clearly holds by independence of 
matrix elements and our zero mean assumption. In the static Hermitian case, one 
needs in addition that E^Zfj = 0. The time-dependent analog reads ¥,Zi2{s)Zi2{t) = 
for any s,t G R, which is one of our assumptions. This allows the exact same 
reduction to go through in the time-dependent setting. 

The only difference of the multi-set case from the one-set case is that one needs 
to keep track of the alphabets the words are built from: A word corresponding to 
coordinate number p of (2) would have to be built from the alphabet that coincides 
with the set Bp. Equivalently, the corresponding graphs will have their vertices 
labeled by elements of Bp. 

Since all sizes \Bp\ have order L, and \Bi U • • • U Brn\ = 0{L). and also the mo- 
ments of matrix elements at all times are uniformly bounded, the estimate showing 
that all contributions not coming from matchings are negligible (Lemma 2.1.34 in 
[AGZ]) carries over without difficulty. It only remains to compute the covariance. 

For real symmetric Wigner matrices in the one-set case the limits of the variances 
of the coordinates of (2) are given by (2.1.44) in [AGZ]. It reads (with k = kp for 
a p between 1 and m) 



(4) 2ecU+k'Cl+J2 — 



r=3 



ki>0 »=1 



where {Ck}k>i are the Catalan numbers, and we assume Ca = unless a G 
{0,1,2,...}. The Catalan number Ck counts the number of rooted planar trees 
with k edges, and different terms of (4) have the following interpretation (see [AGZ] 
for detailed explanations): 

• The first term comes from two trees with (k — l)/2 edges each that hang from 
a common vertex; the factor fc^ originates from choices of certain starting points 
on each tree united with the common vertex, and the extra 2 is actually EY^. 

• The second term comes from two trees with fc/2 edges each that arc glued 
along one edge. There are fc/2 choices of this edge for each of the trees, there is an 
additional 2 = EZfj — 1, and another addional 2 responsible of the choice of the 
orientation of the gluing. 

• The third term comes from two graphs each of which is a cycle of length r 
with pendant trees hanging off each of the vertices of the cycle; the total number of 
edges in the extra trees being (fc — r)/2 (this must be an integer). As for the first 
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term, there is an extra k"^ = k ■ k coming from the choice of the starting points and 
also an extra 2 for the choice of the gluing orientation along the cycle. 

For each of the three terms the total number of vertices in the resulting graph is 
equal to fc, and if one labels each vertex with a letter from an alphabet of cardinality 
\B\ this would yield a factor of 



\B\m-^)---m-k+i) = \Bt+om 

Normalization by \B\^ yields (4). 

In the general case, in order to evaluate the covariance 

(5) 



^s-l^ 



(Tr(X^;(ip)) -ETr(X^;(t,))) - ETr(X^'^(t,))) 

in the limit, we need to employ the same graph counting, except for the two graphs 
being glued now correspond to different values kp and kq of fc, and their vertices 
are marked by letters of different alphabets Bp and Bq. 

• The first term gives 2kpkqC kp-i C k,-i for the graph counting, and an extra 

2 2 

IBpHBql . i\Bp\ - l){\Bp\ - 2) . . . {\Bp\ - ^) • {\Bq\- l){\Bq\ - 2) . . . {\Bq\ - ^) 

for the vertex labeling (the factor \Bpn Bg\ comes from the only common vertex). 

Moreover, W^Y^ is replaced by ¥yi{tp)Yi{tq) = c{tp,tq). Normalized by L~ "2 " 
this yields 



k-n 1 kn — 1 



'ZkphqC kp — 1 C kq — 1 C(tp^ tq^bpqbp bq 
2 2 

kp ^ kq 

• The second term has kpkqCkp Ck, from the graph counting and c^„hp hq 

~2~ ~2~ 

from the label counting. In addition, ¥.Zi2'^ — 1 is replaced by EZf2(ip).^i2(ig) — 1 
2c^{tp,tq). The total contribution is 

i2 it — 1 ^2.-1 
kpkqC kp C kq \^c(tp^tq^bpq ) bp" bq" 



For the third term in the same way we obtain 



E2kpkq 
T 

r=3 



Si>0 

y2 "YLi—i Si— hp 



I 



(^c{tp,tq)bpq^ bp 



I 



where c'^{tp,tq) appeared as {EZi2{tp)Zi2{tq)y , which in its turn came from the 
edges of the r-cycle. 

Thus, the asymptotic value of the covariance (5) is 

2kpkqC kp—2_C kq—2_ic(tp^tq^bpq)bp ^ bq ^ -^kpkqC kp^C kq^{c(tp^tq^bpq) b^ bq 



E2kpkq 
r 

r=3 



E n^^ 



Si>0 »=1 
\2 J2i=l Si=kp-r 



E n^*- 

ti>0 »=i 
\2j:i=iU=kq-r I 



{c{tp, tq)bpq) bp = bq 



We now use the fact that for any 5 = 0, 1, 2, . . . 



r \ r 



Si>0 ' 



V S j2S + r 



see (5.70) in [GKP]. This allows us to rewrite the asymptotic covariance in terms 
of binomial coefficients: 



{kp - 1)/2J V(fc, - l)/2 

kn 



k-,,-1 k„-l 



+ 4 

OO 

r=3 



A;, 



(fcp - r)/2j \{kg - r)/2 



j {c{tp,tq)bpg) bp ' 6q 

{c{tp,tq)bpqYb^ 6g~ 



{c(tp,tq)bpq) bp ^ 6, 



Using the binomial theorem, we can write this expression as a double contour 
integral 
(6) 



II 



z + 



w + 



bqY' C{tp,tq)b. 



w 



PI 



dzdw 



consti — \z\<\w\—const2 



C{tp,tq)bp 



Consider the right-hand side of (3) and assume that \z\'^ = bp < bq = 
Observe that 



2 In 



c(tp,tq)bpq - ZW 



C{tp,tq)bpq - ZW 



■z — w 



= - In 



■ z — w] +ln 



= -2 In 
bp 

+ ln 



This allows us to rewrite the right-hand side of (3) as a double contour integral 
over complete circles in the form 



C{tp,tq)hp 



C{tp, tqjbpq _ . 

— i-—^ — — z — w \ — In 



C(tp,tq)bpq_ _\ 

— H-^ — — z — w ] . 



|^|2 = 6p \w\^=bg 

Recalling that /3 = 1 and noting that 



1 1 / C.{tp, tq)bp 



-z — w 



dx{z) dx{w) 
dz dw 



dzdw. 



kp{x{z))''^ 



_i dx{z) _ djxjz))'''' , , sska-i Mw) _ d{x{w))'''" 

— ! Kq{X[W)) — 



dz 



dw 



dw 



we integrate by parts in z and w and recover (6). The proof for for hp = bq is 
obtained by continuity of both sides, and to see that the needed identity holds for 
bp > bq it suffices to observe that both sides are symmetric in p and q. 

The argument in the case of Hermitian Wigner matrices is exactly the same, 
except in the combinatorial part for the first term the factor 2 is missing due to 
the change in EYi{s)Yi{t), in the second term 2 is missing due to the change in 
E|^i2(s)P|-^i2(i)P, and in the third term 2 is missing because there is no choice in 
the orientation of two r-cyclcs that arc being glued together. □ 

Proof of Proposition 1. We need to show that for any complex numbers {uk}^^i 



y2 ^'^"^ / / 



fk{z,s)fi{w,t)C{z,s;w,t)dzdzdsdwdwdt > 0. 



We can approximate the integration over the three-dimensional domains by finite 
sums of one-dimensional integrals over semi-circles of the form \z\ = const, s = 
const. On each semi-circle we further uniformly approximate the (continuous) 
integrand by a polynomial in ^{z). Finally, for the polynomials the nonnegativity 
follows from Theorem 2'. □ 

Chebyshev polynomials. One way to describe the limiting covariance structure 
in the one-matrix static case is to show that traces of the Chebyshev polynomials 
of the matrix are asymptotically independent, see [J]. A similar effect takes place 
for time-dependent submatrices as well. 

For n = 0, 1, 2, . . . let T„(a;) be the nth degree Chebyshev polynomial of the first 
kind: 

Tn{x) = cos(n arccosx), T„(cos(x)) = cos(na;). 
For any a > 0, let T^{x) = T„(|) be the rescaled version of T„. 
Proposition 3. In the assumptions of Theorem 2", for any p,q = 1, . . . , m 



Urn E 

L—i-oo 



X (TV(r,^/^(XB,(i,)))) -ETr(r,^y^(x^,(i,))) 



2/3 \ 



Proof. Using (6) and assuming bp < bq we obtain that the needed limit equals 

2 f f „ , , . , c{tp,tq)bpq dzdw 



j I Tfe^(cos(arg(^;))Tfc^(cos(arg(w)) 



bp=\z\<\w\=b. 
1 



2/3(27ri)2 



bp = \z\<\w\=bq 




C{tp,tq)h 



pg 



dzdw 



' C(tp,tq)bp 



— w) 



2 • 
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Writing — w) ^ as a series in z/w we arrive at the result. Continuity and 

symmetry of both sides of the Hmiting relation removes the assumption bp < bq. □ 

Note that in the Gaussian specialization (when c{s,t) = exp(— |s — t\)) and 
for a single size L time-dependent Wigner matrix (i.e. bp = bq = bpq = 1), the 
centralized traces of Chebyshev polynomials of this matrix evolve as independent 
Ornstein-Uhlenbeck processes with speeds equal to the degrees of the polynomials. 
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